Surface electronic scattering in d-wave superconductors. 



A. A. Golubov 1 , M. Yu. Kupriyanov 2 
1) Department of Applied Physics, University of Twente, 
P.O.Box 217 
7500 AE Enschede, The Netherlands 
2) Nuclear Physics Institute, Moscow State University, 119899 GSP 
Moscow, Russia 

Theoretical model of a rough surface in a d-wave superconductor is studied for the general case 
of arbitrary strength of electron scattering by an impurity layer covering the surface. Boundary 
conditions for quasiclassical Eilenberger equations are derived at the interface between the impurity 
layer and the d-wave superconductor. The model is applied to the selfconsisrent calculation of the 
surface density of states and the critical current in d-wave Josephson junctions. 
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Understanding the transport properties of surfaces and interfaces in high temperature superconductors (HTS) is 
important for correct interpretation of the experimental data and for applications of these materials. Up to now 
surface peculiarities in d-wave superconductors were extensively discussed in the framework of theoretical models 
based on quasiparticle reflection from specular interfacesata. The influence of surface roughness in this accroach 
was studied by the Ovchinnikov methocfl, namely, by introducing a thin impurity layer covering the surfacd3i3. The 
degree of roughness is then measured by the ratio of the layer thickness d to the quasiparticle mean free path in the 
layer I. This approach was used to study the smearing of Andreev surface bound states by weak disorder (d <C £). 

The regime of strong surface roughness, d ^> i, was discussed irfl. In this case the parameter-free boundary 
conditions were derived at-the interface between the superconductor and the disordered layer. Similar conditions 
were also derived earlier inS at rough walls in superfluid 3 He — B. As was shown iro, an isotropic order parameter 
is nucleated in the disordered layer even in the absence of the bulk subdominant pairing interaction in the s-wave 
channel. This leads to the formation of the gapless s-state at the interface. This phenomenon was called an anomalous 
proximity effect between a d-wave superconductor and a disordered layer. 

In the present paper we consider the anomalous proximity effect between a d-wave superconductor and a surface 
impurity layer of thickness d -C min {V£o^j £o} m the general case of an arbitrary electronic mean free path £. 

We consider the interface oriented normally to the crystallographic ab plane and assume elastic Born scattering. 
The normal and anomalous quasiclassical propagators, g and / , depend on the coordinate x in the direction of the 
surface normal and on-,the angle between the Sjtujface normal and quasiparticle trajectory and obey the quasiclassical 
Eilenberger equationsa. Using the substitutionst3 

f= 2a + 2b = l~ab 
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we rewrite these equations in terms of the scalar differential equations of the Ricatti type 
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b(x) = -A(x)-±-(f+), (2) 



VCOS0 



da(x) 
dx 



2<j n + -(g) + a(x)(A(x) + ^-(f)) 

T ZT 



a(x) = A(x) + — </) . 
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The pair potential A should be found from the selfconsistency equation 



Aln^ + 2^r£(--<A(M , )/>)=0. 
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Here uj n — irT(2n + 1) are the Matsubara frequencies, v is the Fermi velocity, r = Ijv. Assuming a cylindrical Fermi 
surface, angle averages are defined as (...) = (1/2tt) J {...)d9. 

We consider pure d-wave interaction in the bulk with the coupling constantllJ 



A(0, 0') = A d (0, 9') = 2A cos(2(6> - a)) cos(2(6>' - a)), 
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where a is the misorientation angle between the crystallographic a axis and the surface normal. In this case the bulk 
anomalous propagator has the d-wave symmetry 



,/ ^ V2A 00 cos(2(6» - a)) , , 

b(+oo, 6) = a(+oo, 7T + 6) = ; °° v v " 5 

lu + v /w 2 + 2A^ cos 2 (2(#-a)) 

As will be shown below, this symmetry is violated near the interface. 

Equations (§)-(§) must be supplemented with the boundary conditions connecting functions b(0,9) and a(0,6) 
which describe respectively the electrons incident and reflected from the interface. This can be done by matching the 
solutions in the impurity layer (— d < x < 0) and in the clean d-wave region (0 < x). 

Since the thickness of the impurity layer d -C £ e ff — min {\/£o£, £o} ; w ^ can neglect the terms proportional to 
lu and A in (j^) and (||) and consider (/) and (g) as,-spatially-independent constants, which must be determined 
selfconsistently. Making use of the boundary conditions at the totally reflecting free specular interface (x = —d) 

b(-d,-6) = a(-d,d), (6) 

we find the solution of the equations (0), (|3|) at — d < x < in the form 
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From (Q), (||) after a simple algebra we arrived at the boundary condition at the interface between the clean d— wave 
superconductor and disordered layer 

,7fn 0)=h(0 0) (l-Gtanh(fccQ) tanh(kd) 

V ' ; 1 ' j (Fb(0, -9) + G) tanh(fcrf) + 1 (Fb(0, -6) + G) tanh(fcrf) + 1 ' 1 ' 

Equation (||) describes the fact that the amplitude of anomalous Green's function along outgoing trajectories, a(0, 9), 
consists of two parts. The first, specular part, is proportional to 6(0, —6) and describes the correlation between the 
incoming —9 and reflected 9 trajectories. The second, diffusive part, is proportional to F and describes an average 
contribution of all the incident trajectories to the outgoing one in the 9— direction. 

In the limit of weak disorder (d <C £) the boundary condition (||) reduces to the specular one (|J). At finite d there 
is a cone of angles arccos(d/€) < 9 < ir/2 in which the scattering from the interface is rather diffusive than specular. 
The larger is d the smaller is the correlation between incoming and outgoing trajectories. n 

In the limit of strong disorder (d ^> £) the solutions (f?j)- (||) transform to those of Usadel equations obtained inEl 
and it follows from (||) that 

a(0) = i +{9M y m 

This is exactly the condition used previously i for the analysis of strongly disordered d-wave interfaces. According 
to this condition the incoming and outgoing trajectories at the rough interface with strong disorder are completely 
uncorrelated as a result of the scattering from the impurity layer. 

It is straightforward to find selfconsistent solutions of the set of equations (||)-(f|), @- First, we have numerically 
integrated the equation (|^) for bjx) with the initial condition (^|) along the trajectory from x = +00 (bulk) to x = 
(interface). Then the equation ([|) for a(x) was integrated from x = to x — +00 with the initial condition (^|). The 
angle averages (/(0)) and (g(0)) and the pair potential A(x) were calculated iteratively. 

We have applied this model for calculation of the low-temperature conductance of NID tunnel junction, which is 
expressed through the surface density of states N(0,9;e) 
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J_J /2 D(0) cos OdO 

where Rm is the normal state junction resistance and D(9) is the barrier transmission coefficient. The angle-resolved 
surface density of states is given by N(0, 9; e) = Re(g(0, 8;uj = —is))- 

Fig.l shows the results of calculations for a = 20° and various values of the surface scattering parameter d/£. We 
have chosen <5-shaped potential barrier with the angular dependence of transmission D(9) = D(0) cos 2 0. Zero-bias 
anomaly (ZBA), which is the manifestation of the sign reversal of the d-wave pair potential, is strongly.-broadened 
with increase of d/£. In the limit of weak disorder d/£ <C 1 the results exactly correspond to those of Rem. 

With the increase oid/i ZBA is smeared out completely, while the signatures of finite-energy Andreev bound states 
are still present. The latter are due to quasiparticles trapped in the surface region with the reduced pair potential 
A(a;). As is seen from Fig.l, the angle-averaged density of states in the disordered layer is gapless and has a number 
of peaks at the energy below the bulk pair potential. Note that for large values of djl there are several (weak) 
peaks corresponding to more than_pne bound state, in contrast to the case of specular interface when only a single 
bound state exists at finite energju. This is because the average trajectory length of trapped quasiparticles, which 
determines the number of resonances, is larger than £o in the diffusive case. Mathematically this fact is described by 
the boundary condition Eq.|[ 

Next we apply the model to the study the influence of surface roughness on the Josephson supercurrent .in.a .tunnel 
junction based on d-wave superconductors (DID junction). The critical current is given by the expressiontailj 

nT ^E u >o/-^ 2 [/i.W/2-W- ha{6)ha{6)]D{6)cos6d6 

eR N I c = —j- . (12) 

Jli /2 D(9) cos 9d9 

Here fi^s(a) — (/1.2 ± ft 2)!^ ancl /i,2 are the surface Green's functions of the left (right) electrode. 

The critical current was calculated numerically using the solutions for f 12 and assuming D{9) = D(0) cos 2 9. We 
consider two cases: (1) symmetric junction with the misorientation angles 0.1 = 0.2 = 20°, (2) antisymmetric (mirror) 
junction wifli a\ = —a 2 = 20°. The results of calculations are shown in Fig. 2. In complete agreement with the results 
of Ref.trta the mirror junction exhibits a nonmonotonous / c (T)-dependence due to the anomaly /1,2a ~ at low 
frequencies related to the midgap surface bound state. Like ZBA in conductance, this anomaly is smeared out with 
the increase of the surface scattering parameter d/£. 

In the strong scattering regime d/^> 1 the I c Rn product of the DID junction becomes small and saturates at the 
level el^ is Rn /2nT c ~ 1CP 3 (the corresponding number for a conventional SIS junction is e/f is Rn/2ttT c ~ 0.44). 
The anomalous contribution f\ a f2a vanishes in this case and I c for the symmetric and mirror junctions coincide. 
However, as shown in inset in Fig. 2, I C (T) becomes nonmonotonous with an increase of djl. The reason is the 
destructive interference within the impurity layer: the phases of Eilenberger functions on the incoming trajectories 
alternate due to the d-wave angular structure. As a result, the angle average (/1.2) vanishes at oj n <C ttT c , reaches 
a maximum at uj n ~ irT c and falls down as lu~ 2 at uj n 3> ttT c . It is the property /1,2s - * at small uj n which is the 
reason of decreasing I C (T) at low T. Variation of the misorientation angle a for djl 3> 1 does not change the I C (T) 
behavior, while I c Rn product at T = scales as cos 2 (2a). 

In conclusion, we have derived the boundary conditions for the Eilenberger functions at the rough surface of a 
d-wave superconductor and applied them to the study of the crossover from specular to diffusive surface scattering. 
The low-temperature conductance of NID tunnel junction and the Josephson supercurrent in DID tunnel junction are 
calculated for arbitrary degree of surface roughness. It is shown that the width of the anomalies in the conductance 
and in the critical current is controlled by a single scattering parameter d/i. 
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FIG. 2. Influence of surface roughness on the temperature 
dependence of the critical current in the DID tunnel junc- 
tion. Dashed lines: symmetrical junction cti = ct2 = 20°. 
Dotted lines: mirror junction an = — Q2 = 20°. Inset: non- 
monotonous behavior of I c for d/£ 3> 1. 



